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We examine the problem of simulating lattice gauge theories on a universal quantum computer. 
The basic strategy of our approach is to transcribe lattice gauge theories in the Hamiltonian for- 
mulation into a Hamiltonian involving only Pauli spin operators such that the simulation can be 
performed on a quantum computer using only one and two qubit manipulations. We examine three 
models, the U(l), SU(2), and SU(3) lattice gauge theories which are transcribed into a spin Hamil- 
tonian up to a cutoff in the Hilbert space of the gauge fields on the lattice. The number of qubits 
required for storing a particular state is found to have a linear dependence with the total number 
of lattice sites. The number of qubit operations required for performing the time evolution corre- 
sponding to the Hamiltonian is found to be between a linear to quadratic function of the number 
of lattice sites, depending on the arrangement of qubits in the quantum computer. We remark that 
our results may also be easily generalized to higher SU(iV) gauge theories. 

PACS numbers: 03.67.Lx,11.15.Ha,12.38.Gc 
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I. INTRODUCTION 

The efficient simulation of quantum many-body prob- 
lems is today of vital importance to many fields 
of physics, stretching from condensed matter physics, 
atomic physics, molecular and optical physics, nuclear 
physics, quantum chemistry, and high-energy physics. 
Despite advances in various numerical techniques, a truly 
reliable, accurate, and efficient method for extracting 
physical quantities for such problems does not exist in 
general. For example, quantum Monte Carlo methods 
suffer from the "minus-sign problem," and the more re- 
cent Density Matrix Renormalization Group (DMRG) 
method possesses difficulties in dimensions greater than 
one 0. The impact of a method that is free from such 
difficulties would undoubtedly offer new insights into the 
physics of quantum many-body systems. 

Quantum computers are hoped to offer an alternative 
to such methods. The exponential parallelism that is 
employed by quantum computers is believed to overcome 
the problems associated with the exponential explosion 
of Hilbert space size inherent to quantum many-body 
problems. It has been shown that an efficient simula- 
tion is indeed possible for a large class of spin, fermionic, 
and bosonic Hamiltonians yi-Ll- For example, fermionic 
Hamiltonians are implemented via an algebraic mapping 
of operators via the Jordan- Wigner transformation. This 
allows fermion operators appearing in the Hamiltonian 
to be written entirely in terms of spin operators. As 
discussed in Ref. @], this approach is free of "minus- 
sign problem" that is present in Monte Carlo methods. 
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The simulation is then performed by evolving the system 
according to the time evolution operator by a series of 
gates which evolve the system according to the system 
Hamiltonian. A phase estimation algorithm using the 
quantum Fourier transform or a similar algorithm is then 
performed to calculate the eigenenergies of a given sys- 
tem We assume implicitly here the standard model 
of quantum computation with two-level spin qubits and 
unitary operations. For bosonic operators, there is no ex- 
act algebraic mapping of operators into spin operators, as 
can be simply understood 5] by the disparate sizes of the 
Hilbert spaces - a single boson has an infinite dimensional 
Hilbert space, whereas a spin has only two dimensions. 
However, it is still possible to map the bosonic Hamil- 
tonian onto a spin Hamiltonian via a transcription of 
the operator mappings 0. Clearly due to the mismatch 
in the system dimensions of the bosons and spins, there 
must be some truncation in the Hilbert space. However, 
this truncation can have little or no effect in terms of the 
Hamiltonian if performed in the correct way. For exam- 
ple, Hamiltonians that conserve particle number have a 
maximum of the total number of particles on a given site, 
and therefore has an effective cutoff equal to the number 
of particles. 

One problem that has received little attention is the 
problem of the quantum simulation of lattice gauge the- 
ories 0, Q . Arguably this is the most computationally 
intensive quantum many-body problem of all, due to the 
large numbers of degrees of freedom per site and the ne- 
cessity of simulating in three spatial dimensions. Since 
the first Monte Carlo simulations over 25 years ago 
the lattice gauge theory program has steadily advanced 
with increases in computer speed and improvements in 
simulation techniques. For the first time in the last few 
years it has become possible to perform a realistic simu- 
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lation full QCD (i.e. with fermions) without quenching 
;9j. However, should a method become available that 
could give accurate results while being computationally 
less expensive, this could offer a valuable alternative to 
the methods being employed today. 

In this paper we study the problem of whether it is 
possible to implement lattice gauge theories on a uni- 
versal quantum computer efficiently. Before commencing 
we must make several choices in how we formulate the 
problem. First, we assume that a quantum computer 
with individually addressable qubits and controllable lo- 
cal (nearest neighbor) two-qubit interactions is available. 
In particular, we assume that our quantum computer 
consists of two-level qubits, so that it is operated on by 
Pauli spin operators. This is simply for convenience as 
most of the literature is formulated in this way. Sec- 
ond, we choose to concentrate on pure gauge field the- 
ories alone, as an extension to include fermions should 
involve only a small modification to the results presented 
here. For example, fermion fields may be implemented 
in a similar way to that described in Ref. |2j. The third 
regards the formulation of the lattice gauge theory. To 
take advantage of the quantum parallelism of the quan- 
tum computer we choose a Hamiltonian formulation of 
lattice gauge theory, in contrast to the Euclidean for- 
mulation that is more commonplace in present-day sim- 
ulations. The basic strategy is then to rewrite the op- 
erators appearing in the lattice gauge Hamiltonian into 
terms of spin operators. As is the case for bosons, no 
algebraic mapping of operators to spin operators is avail- 
able for gauge operators due to different Hilbert space 
sizes. Analogously to boson operators, we will show that 
it is possible to preserve the operator mappings of the 
gauge operators. After this is done, one can implement 
the methods described in Refs. H 113 EH to perform 
the evolution. The Hamiltonian formulation has also the 
advantage that the lattice dimension for a simulation of 
real-world QCD (for example) is three rather than four 
for the Euclidean formulation, which reduces the space 
and time requirements for the calculation. 

We consider three models in this paper. First we con- 
sider compact U(l) lattice gauge theory. Although the 
transcription is particularly simple in this case, this will 
serve to introduce the basic strategy that is employed. 
We then move onto the next most complicated case of 
SU(2) lattice gauge theory. This contains all the ele- 
ments that are required in order to transcribe any Hamil- 
tonian for SU(A). We finally consider the transcription 
of the SU(3) Hamiltonian, which is the closest relation 
to QCD. Although it is crucial that the time evolution 
operator can be implemented efficiently, it is also impor- 
tant that other steps in the simulation can be carried out 
efficiently. To show that this is possible, in Sec. IIVI we 
discuss methods for initializing the qubit states, suitable 
for extracting quantities such as the eigenenergies and 
expectation values of the low-lying eigenstates. This is 
again carried out for the U(l), SU(2), and SU(3) lattice 
gauge theories. 



This paper is organized as follows. In Sec. [H]we give a 
brief review of the simulation of spin Hamiltonians (Sec. 
Ill Aj) and lattice gauge theory (Sec. IIIB|) . In Sec. HTH 
we transcribe the compact U(l) lattice gauge theory (Sec. 
IIII Af) . SU(2) lattice gauge theory (Sec. IIIIBfl . and finally 
the SU(3) lattice gauge theory (Sec. 1111 Cfl Hamiltonians. 
Sec. IIVI discusses the initialization of qubits and Sec. El 
discusses how observables may be extracted. Sec. IVII 
gives a summary of our findings and conclusions. 

II. PRELIMINARIES 
A. Simulating spin Hamiltonians 

We now give a brief review of the method of simu- 
lating spin Hamiltonians following Refs. 0, The 
various methods of extracting eigenstates, eigenenergies, 
and expectation values for a given Hamiltonian in these 
references rely upon the ability of performing the time 
evolution corresponding to U(t) — ex-p(—iHt) where H 
is the Hamiltonian under investigation 0, 0, E3, 0, 0] ■ 
We cannot directly perform the evolution U{t) since we 
assume only one and two-qubit manipulations are avail- 
able on our quantum computer. We therefore decompose 
the evolution operator using the Trotter formula jjjj]: 

e ~iHt _ ( e -iH 1 t/m e -iH 2 t/m e ~iH NH t/m\ m 

m — >oc V / 

(1) 

where the Hamiltonian is a sum of Nu terms H = 
Yli=i Hi ■ Assuming that the Hamiltonian under consid- 
eration here is entirely composed of Pauli spin operators, 
each term on the RHS of Q will consist of the exponen- 
tial of a product of a number of spin operators 

Ui(t)=eM-iTil[<r?), (2) 

j 

where 7, is a constant involving coupling constants and 
the time of evolution, j is an operator label (e.g. labeling 
the site number etc.), and aj = x, y, z. Operators such as 
(0) may be decomposed entirely into one and two-qubit 
terms. Denote single qubit manipulations by 

R j {e,n) = e i(e ^- n+s \ (3) 

where S is a global phase shift. Denote two-qubit manip- 
ulations by Rjk(uj). For example, a two-qubit manipula- 
tion that is available may be 

R jk {u) = e-^°i. (4) 

The precise form of the two-qubit interaction is not im- 
portant, so long that one is present. Terms that involve 
only one spin operator in the exponent in (0) may be 
directly constructed from Terms involving a prod- 
uct of two or more spin operators in the exponent of @ 
can be written by an appropriate combination of one and 
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two-qubit operations. For example, say we need to con- 
struct the operator e~ l(Tl a ? . This may be decomposed 
into the product 

e -iata\ = R l^/^ x ) R \^/± iV ) Rl2 (i) Rl ^/A,y)R 2 {n/^x) 

where we assumed the two qubit-manipulation and 

used the identity We M U = e u MU , where U is an uni- 
tary operator and M is an arbitrary operator. Longer 
products of spin operators may be constructed by suc- 
cessive products of two-qubit operators Q • 

The simulation then proceeds in the following way. 
The initial state of the qubits is first prepared, the pre- 
cise form of which depends on the quantity that is being 
calculated. This initial state is then evolved according 
to the particular algorithm to be performed, which im- 
plements the time evolution operator U(t) — e~ lHt . The 
time evolution is performed according to the method de- 
scribed above. Finally a measurement is performed on 
the system to extract the desired results. For example, 
in order to measure the eigenvalue spectrum, the qubits 
are prepared in a state with a non-zero overlap with the 
eigenstate of interest. This is typically the ground or low- 
lying eigenstates of the system. A suitable initial state 
is for example the mean-field ground state of the Hamil- 
tonian. This initial state is then evolved forwards using 
the time evolution operator U(t), using a phase estima- 
tion algorithm [l^. The estimated phases then recover 
the eigenspectrum of the Hamiltonian. Further details 
of the extraction of observables may be found in Refs. 
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FIG. 1: An elementary plaquette. Each link is labeled by 
the coordinate r of one of the sites that it connects, and the 
direction of the link [i, v. 



and 



ifocPl T l^ ( 7 ) 

where 5 a @ is the Kronecker delta and f af3 "> are the struc- 
ture constants of the group. The following commutation 
relations then follow 



k=l 
N 



[E c 



I] 



= if 



k=l 



(8) 



(9) 
(10) 



where the indices i,j on Uij refer to the matrix elements 
of the N x N matrix. We denote all quantities which 
contain the operators E a or A a in bold font. The SU(iV) 
pure gauge Hamiltonian is written 



B. Lattice gauge theory 



Let us define SU(AQ lattice gauge theory in the Hamil- where 



m 



tonian formulation |13( (an excellent review of lattice 
gauge theory in the Hamiltonian formulation may be 
found in Ref. 0])- Assume a <i-dimensional lattice. 
Label the sites of the lattice by a coordinate r. Each 
site has 2d links emanating from it. Label each of these 
links according to the site it is attached to, and the di- 
rection that the link points (see Fig. ^for the labeling 
conventions). On each link of the lattice define an TV 2 — 1 
component operator A", with canonically conjugate op- 
erators E a . These are the gauge field operators of the 
lattice gauge theory. Define an N x N unitary matrix 
involving the gauge fields 



U(r,fi) = exp(ziga 



M°(r,A0) 



(5) 



where r Q (a — 1, ... , N 2 — 1) are the fundamental gener- 
ators of SU(iV), g is the coupling constant, and a is the 
lattice spacing. The generators of SU(iV) obey 



1 ™/3 
2 U 



(6) 



(2N-(Z( P )+Z\p)))} (11 



pG {plaqucttos} 



Zip) = Tr [U(r, fi)U(r + fi, v)U\r + v, ^)U\r, v)] 

(12) 

and x = Jr. The trace in the Z(p) operator is taken after 
matrix multiplication of the four U matrices. The first 
summation in the Hamiltonian (|llfl sums over all links in 
the lattice, and the summation involving Z(p) operator 
sums over all plaquettes in the lattice. The arrangement 
of links in a plaquette operator is shown in Fig. ^ The 
overall constant of g 2 /2a and the constant term 2N in 
the plaquette summation merely rescales and shifts the 
energies respectively and will be omitted for the rest of 
this paper. 



III. FORMULATING LATTICE GAUGE 
THEORIES VIA SPIN OPERATORS 

A. Hamiltonian for U(l) lattice gauge theory 

We first start with one of the simplest lattice gauge the- 
ories: compact U(l) lattice gauge theory (compact quan- 
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turn electrodynamics with no fermions). The transcrip- 
tion of the operators is particularly simple in this case 
and will therefore serve to introduce the basic method 
that will be employed for the SU(2) and SU(3) cases. 
The U(l) Hamiltonian is pi 



r.//, 



pG{plaquettes} 



Z{p) + Z\p)) 

(13) 

where 

Z(p) = e ie ( r >^) e ie ( r +A I ^) e - ,;£ '(''+ ! ^) e - i£ '( , % 1 ') (14) 

where we have defined 9 = gaA according to standard 
notation. The generator of U(l) is a scalar and therefore 
all the U = e commute with each other. The U are 
lxl matrices. The E obey the commutation relations 



[E,U] = U 
[E,U*] = -UK 



(15) 
(16) 



The Hilbert space for the Hamiltonian may be generated 
by successive applications of Z(p) + Z<(p) on the state 
|0) which is defined to be an eigenstate of the E operator 
satisfying 

£(r,/i)|0> = (17) 

for all r and fi 0] . 

Let us first consider the Hilbert space of a single link. 
A link where U has been applied I times is an eigenstate 
of E 2 with eigenvalue I 2 : 

E 2 (U) l \0) =l 2 (U) l \0). (18) 

We can therefore identify the orthogonal states 

\l) = (U) l \0) (19) 

which satisfy 

E 2 \l) = l 2 \l) (20) 

U\l)=\l + l) rf\l) = \l-l). (21) 

The plaquette operator then shifts the \l) eigenstates on 
links arranged on all squares of the lattice (see Fig. [2j. 

This space of states may be easily rewritten in the 
spin language in the following way. On a given link, a 
quantum register of D qubits keeps track of the \l) state 
that the link is in. Denote the state that the register is in 
by |i) ta g, where the "reg" serves as a reminder that this 
is a state of the register. Clearly, since I = {— oo, . . . , oo} 
and D qubits can only keep track of a maximum of 2 D 
values I, some truncation is involved. This truncation 
is performed such that I = {— l max , . . . , Z max }. D qubits 
will then give an £ max of (2 D — l)/2. The impact of this 
truncation can be handled by analyzing the convergence 
of observables (e.g. the energy spectrum) with l max |3lj . 
Now define operators that act on these registers 

E 2 \l)^ = l 2 \l)^ s (22) 

L±\l) tes - |Z±l) teg (23) 

L + \L ax ) rcs = (24) 

L-l-hJ)^ = (25) 



/=-! 



/=-l 



1=1 



1=1 



FIG. 2: A single plaquette excitation for U(l) lattice gauge 
theory. All unlabeled links have / = 0. 



A simple implementation of these operators is given in 
Appendix ^ It is then a simple matter to rewrite the 
U(l) Hamiltonian as 

H vm =^E*(r,n)-x £ (Z(p) + ZHp)) (26) 

r,/l p£ {plaquettes} 

where 

Z(p) = L + (r, f i)L + {r + (i,v)L-(r + is,Li)L-(r,v), 
Z\p) = L-(r,n)L-{r + (i,v)L + (r + v,[i)L+(r,v), 

A similar Hamiltonian was given in Refs. 0,^3- 

For a total of M lattice sites, this requires a total of 
~ MdD qubits and is therefore linear in space require- 
ments. Let us estimate the dependence of the total num- 
ber of operations required for a particular time evolution 
U(t) with respect to the number of lattice sites M. By 
a judicious arrangement of qubits, it is clear that the M 
dependence can be made linear. In Fig. |3| the qubits 
associated with an \l) register for a particular link are 
grouped spatially in the same location, with the same 
units repeating in the same configuration as the links 
on the original lattice. As each term in the Hamilto- 
nian l|26|) only acts between nearest neighbors, there is 
no dependence on M for each term. The total number 
of operations is therefore proportional to the number of 
terms in the Hamiltonian, which is cx Aid. A less ideal 
situation will have a qubit arrangement which involves 
manipulations between qubits which are spatially sepa- 
rated by distances of the order of the number of qubits 
required for the simulation. In the worst case, every term 
in the Hamiltonian will involve qubit operations between 
very distant qubits. The number of operations, i.e. the 
sequence of swap operations, required for executing such 
a long-range term is of the order of the total number of 
qubits ~ MdD. Therefore, the total number of opera- 
tions will grow oc M 2 d 2 D. This is quadratic with the 
lattice size. The best case and the worst case scenarios 
presented here are both probably rather unrealistic; we 
expect that the true dependence will fall somewhere in 
between. 
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[Z(p)f |0>= 




FIG. 3: An optimal arrangement of qubits for lattice gauge 
theories. The qubits associated with a particular link are 
shown in boxes. Registers keeping track of the state of a 
particular link are spatially arranged in the same way as the 
links of the original lattice. 



FIG. 4: A double plaquette excitation, 
have E a (r,n)\0) = 0. 



All unlabeled links 



using the group theoretic identity [l8L 



W mLmR {ct>)W nLnR {4>) 



E 

j=\3-r 



where 



(J, M L \j, n L ;j', m L ) (J, M R \j, n R ;f, m R )Ui lLMR {(j)), 

(29) 



B. Hamiltonian for SU(2) lattice gauge theory 



We now examine SU(2) lattice gauge theory. The 
Hamiltonian is [Dl 



h sv{2) = E ^)) 2 - 2x E z( v) ( 2? ) 



r,fi 



Qfz {plaqucttcs} 



where Z{p) was defined in \\'2\. A simplification oc- 
curred in the above Hamiltonian due to Z^(p) — Z(p) 
for SU(2). The U operators appearing in the plaquette 
operator Z(p) are now 2x2 matrices. The generators 
r a in (JSJ) are the 2x2 Pauli matrices. The commutation 
relations of the E a with the U obey relations © - (fTOft . 

The Hilbert space of the Hamiltonian is generated 
by successive applications of Z(p) on all possible loca- 
tions on the lattice starting from the state satisfying 
E a (r, /i)|0) = for all r and /i. Consider a state such 
as that shown in Fig. 21 where two plaquette operators 
Z(p) operate on the same plaquette. Concentrating on 
one of the four links, this creates terms such as 



U mLmn U nLna |0), 



(28) 



where there is no matrix multiplication between the C/'s 
as all indices have been contracted in the trace of each 
Z(p) operator. We distinguish between the left and right 
matrix element labels using the subscripts L and R on 
the matrix element labels. Such terms may be simplified 



exp(i^TjV a ) 



(30) 



and t" is the generator for the irreducible representa- 
tion of SU(2) corresponding to angular momentum j (t? 
is a (2j + 1) x (2j + 1) matrix). The (J, M\j, n;j', m) 
denote Clebsch-Gordan coefficients for SU(2). We may 
now evaluate l|28|l. Let us generalize this a little by tak- 
ing U 1lLnR — > n . Since the U operators appearing 
in Z{p) are all in the j — 1/2 representation, this corre- 
sponds to j' = 1/2 in and thus there are only two 
terms in the summation: J = j ±1/2 (unless j = 0). It 
is also convenient to define 201 



\j,n L ,n R ) = ^2j + lUi LnR \0), 



(31) 



where the factor of + 1 necessary so that the states 
are correctly normalized: 



(j',m L ,m R \j,n L ,n R ) = 6 3 j>Sm L n L S mR n R - 
We then have 



(32) 



U^ mrt \j,n L ,n R ) 
i+i/2 



E 



J=\j-l/2\ 



2j + l 
2J+ 1 



(J,M L \j,n L ; ±m L )(J,M R \j,n R ; \,m R ) 



x I J, Ml = rriL + til , M R = m R + n R ) 



(33) 



The difference to the analogous result for the U(l) case 
(see Eqn. (|21|l ) is that there are now two terms and 
Clebsch-Gordan coefficients are present. 

Let us also find the effect of operating the first term in 
the Hamiltonian (|27|l on the state \j, riL,n R ). Using the 
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TABLE I: Formulas for Clebsch-Gordan coefficients {J = j + 
Aj,M = rrii + Am\j,rrn; J, Am), where i = L,R is put in 
as appropriate. The Wigner operators corresponding to the 
shifts are also listed for comparison with Refs. |2^. l23jl. 



Coefficient 


Wigner 


Aj 


Am 


Formula 




operator 









Cl2 



('21 



c 22 



1 

1 

1 
1 

f 



f 

1 



f 





1/2 
1/2 
-1/2 
-1/2 



1/2 
-1/2 

1/2 
-1/2 



j + rrij+l 
2j+l 



2j+l 



2j + l 



2j + l 



commutation relation 



3 

n L n R 



a = l m 



we have 



fS a ) 2 |i,ni„nii) = i(j + l)|j,n L ,n fl ). 



(34) 



Our problem is then to rewrite Eqns. (|33|) and H34|) 
into a spin language. A similar problem was examined 
by Bacon, Chuang, and Harrow [21|, where an algorithm 
was presented to perform the transformation 



\j,m 

J+l/2 

E 

J=\j-l/2\ 



(J, M| j, to; i Ato)|J, M = to + Am), 



(35) 



The Clebsch-Gordan coefficients (J, M|j, to; ^, Ato) may 
be calculated using the formulas presented in Tabled A 
completely general way of calcula ting such coefficients for 
SU(AT) is presented in Refs. [H . The circuit of Ba- 
con, Chuang, and Harrow may be rewritten in operator 
form defining the operators 

M ± \m) I 
For the two cases Ato = 



V(Ato = 1/2) 
V(Am = -1/2) 



= |i±|U ( 36 ) 
= |m±|) reg . (37) 

| and Am = —i, we have 

M + (J + c n + J"C2i) (38) 
M-{J+c 12 + J-c 22 ), (39) 



where the c mn are diagonal operator versions of the c mn 
coefficients that extract the coefficients from |j, m) rog us- 
ing the formulas presented in Table [I] An example of 



how such an operator may be constructed is shown in 
Appendix [S] One may verify that this is equivalent to 
the transformation (|35|l by applying the V^Ato) opera- 
tors to the state \j, m) Iog . 

It is now clear how to write an operation correspond- 
ing to (|33|l . For each link, we require three registers 
corresponding to the state \j, rriL,m R ). Denote this by 
\j, tol, TOfl) rog . As for the U(l) case, a maximum cutoff 
must be imposed on these registers to keep the number 
of qubits finite. For j = {0, |, . . . , j max }, the m L and m R 
registers will be in the range {— j max , — i ma x + i) ■ ■ ■ , jmax}- 

We may then write the operator corresponding to i|33|) 

as 



u mn ^ y mn — al m lvl n 

where we defined 

rcg 

and 

M*\j,n L ,n R ) teg -- 
and also 



\j,n L + hn R ) tB 

\j,TlL - 2)"i?)ro 



\j,n L ,n R + ±), 
\j,n L ,n R - 2 ) r 



N k = 



2j+l 
2j+2 

2j+l 
2j 



k = 1 

k = 2 



(40) 

TO = 1 
TO = 2 

(41) 

TO = 1 
TO = 2 

(42) 



(43) 



The diagonal operators Nk recover the normalization fac- 
tors from the registers \j, n^, n R ) r< , g according to l|43|) . 

We may now write the full Hamiltonian for SU(2) using 
these operators. It is first convenient to rewrite (|12|l using 
the identity for SU(2) g3 



rrf _ L]]m-nn 

^ mn I LJ—m—nj 

which gives 

Z(p) = E (-l) mi - m3 U mim2 {r^) 



(44) 



mim2m^m4 

xU m2m3 (r + n,v)U- 

7713—7714 777.4 —mi 

(r, v) 



Let us also define the operator 

E 2 \j,m L ,m R ) ICS = j(j + l)\j,m L ,m R ) T 
We thus have 



(45) 
(46) 



r,/J, p£{plaquettes} 

where 

Z (P) = E {-l) mi - m3 V mim2 (r^) 

mim2m3m4 

x V m2m3 (r + M, v)V- 7713 — 777,4 

(r + is, n)V- 

fii4-mi 

(r, v). 
(48) 
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This is the final result for SU(2). The Hermiticity of this 
Hamiltonian is guaranteed from the property 



where 



vL = (-1) 



(49) 



KM 



exp(i ^ t l u 



0°) 



(52) 



Let us now estimate the number of qubits and the num- 
ber of operations required for this Hamiltonian. Three 
registers are required per lattice link and therefore the 
number of qubits is ~ 3MDd, where D is the number of 
qubits in a single register. This is linear with the num- 
ber of sites M. As for the U(l) case, with a judicious 
arrangement of qubits (see Fig. |3J), the local nature of 
the Hamiltonian may be preserved in the transcribed ver- 
sion of the Hamiltonian. This makes the total number of 
operations cx Md as for the U(l) case. The proportion- 
ality constant will be larger in this case, as each lattice 
link requires a larger number of qubits. There is also 
the added complexity in calculating the Clebsch-Gordan 
coefficients cy, which may be performed in poly(j max ) op- 
erations. Again the worst-case estimate may be obtained 
by assuming that each term in the Trotter expansion is 
a product of operators with a length of the order of the 
number of qubits. We therefore obtain the same result as 
the U(l) case, such that the total number of operations 
is cx M 2 d 2 D. 



C. Hamiltonian for SU(3) lattice gauge theory 

The last case we consider is SU(3) lattice gauge theory, 
which is the theory of principal interest as it is closest to 
QCD. The Hamiltonian is 



and r" are the generators of SU(3) for the representation 
ui. The notation in l|51|l requires some explanation. First 
recall that every state in SU(3) can be labeled by five 
numbers 



\p,q,T,T',Y), 



(53) 



where p and q label the irreducible representation, T and 
T z label the T-spin, and Y labels the hypercharge [24| . 
The label u> in 152fl then corresponds to (p, q) and the 
labels Xl and Xr then label two sets of (T,T Z ,Y). For 
example, the states of the u> = 3 representation (i.e. p = 
1, q = 0) are labeled (in the notation of l |53j l): 

|l,0,i,i,i) «- A = l 



|i,o A 



2 ' 



2' 3/ 



|1,0,0,0,- 



The states of lo — 3* are 



<-> A = 2 
<-> A = 3. 



|0, 1, 2> 2' — l) ^ A — 1 

|0,1, «-> A = 2 

|0,1,0,0,|) <-> A = 3. 



E 

r,/X a—1 {plaquettes} 

(50) 

where Z(p) was defined in (|12|l . The U operators ap- 
pearing in the Z(p) operators are 3x3 matrices. There 
are eight terms in the a summation in JSJ correspond- 
ing to the eight generators of SU(3). The commutation 
relations of the E a with the U obey relations © - (fT0|> . 

The basic strategy is identical to the SU(2) case, al- 
beit with the added complexity of the SU(3) group. Let 
us again consider what kind of states will be generated 
by successive applications of the plaquette terms in lf5U|) 
on the state satisfying E a (r, /x) |0) = for all r and /x. 
Consider again the double plaquette excitation shown in 
Fig. 0] We again obtain terms such as U lLlR U\ L \ R \0), 
which may be simplified using the SU(3) version of the 
identity [lM fl9j: 

(0, T L \u, X L ; J, lL ){^ r K |w, X R ; J, irW^^) 

nr L r R 

(51) 



The (ft, T\oj, A; u', 7) thus denote Clebsch-Gordan cocffi- 
Z(p) + Z\p)) cients for SU(3). We may therefore define 



\w,\l,\r) = \ P ,q,T L ,T[,Y L ,T R ,T^YR) = 

V / MM)^ 1 A K |0>- (54) 

where the factor under the square root is the dimension 
of the irreducible representation (p, q) 

dim(p,q) = (l + p)(l +q )(l+l(p + q)). (55) 

This factor is necessary in order to normalize the states: 



(u> , X L ,X R \u>, X l , X r ) = 8 UU ;'6x L \> L 8x R \> R . (56) 

The U\ L \ R operators appearing in the Hamiltonian (| 5 C)|> 
are all in the fundamental representation of SU(3), i.e. 
the 3 representation (p = 1, q = 0). The U\ x op- 
erators on the other hand are in the 3* representation 
(p = 0, q= 1). Terms such as U^ R U" x \0) may then 
be explicitly written 
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U wm IP. 9> T l Y l^r, T Z R , Y R ) = 

E ^ t/^^y.«',n,22 / .niP,?,T £) 22,li;3 > t £> t£ > y i; ) 

b',«')n=i^-^m=i r «-*Bi v 1 ' 

x (p', q', T R , T Z R \ Y R \p, q, T R , T R , Y R ; 3, t R , t z R , y R ) \p', q', T' L , T£', Y[,T R , T R , Y R ), (57) 



where the (p^q 1 ) summation is over the three terms 
{(P+1 ) «),(P-1,9+1),0,?-1)} and the (t*,tf ,2/i) are 
the corresponding T-spin and hypercharge values corre- 
sponding to the index 7^ (i = L,R). The hypercharge 
and the T z -spin simply add together as 



follows 



Y( = Y l + yi 



(58) 
(59) 



for i = L,R. We also require the rule for operating R 
since the analogous identity to l|44(l does not hold for 
SU(3). This is given by replacing 3^3* in (JSJJ and 
summing over {(j>, q+ 1), (p+ 1, q — 1), (p — 1, q)} in the 
(p 1 , q') summation. 

The first term in (|50|l corresponds to operating the 
Casimir operator for SU(3). We thus have [24J 

V(£; Q ) 2 |p, q) = i (p 2 + q 2 + pq + 3(p + <?)) |p, g>, 



where we have omitted the T-spin and Y labels for 
brevity. 

We may now write the SU(3) Hamiltonian in a spin lan- 
guage. The basic strategy is the same as for SU(2), where 
we use a modification of the circuit of Bacon, Chuang, 
and Harrow |2l|. In Ref |2l| it is suggested that qudit 
registers are used to keep track of the states for SU(iV). 
We follow a slightly different approach, where the same 
qubit registers are used. Therefore all operators acting 
on the registers that we discuss in the following may be 
built from standard SU(2) Pauli matrices. 

For each link associate eight registers according to i|54[l. 
which we denote by \p, q, T L , T[, Y L , T R , T R , Y R ) Ttls ,. As 
in the U(l) and SU(2) cases we must impose a cutoff on 
these registers to keep the number of qubits finite. If the 
highest representation that is stored is then 
we have registers in the range p — {0, f , . . . , p max } and 
q = {0, 1, ... , g max }. The other registers will then be in 
the range 

Tt = {0,i,...,i(p + <?)}, (60) 

t-- = {-§(?+«), -!(?+?) + !,...,!(?+?)}, 

(61) 

Yi = {-§(<? + 2p),-±(g + 2p) + i,...,±(p +2?)}, 

(62) 

for i = L, R. Define operators that shift the registers as 



p ± \ P U s = \p±i) Te 



T?\TA 



Q ± k) I o 

rpZ^Z \ r J lZ 



l<?±l)re g 
= \T Z ± i) 

(63) 



YT\Yi)^ = \Y i ±^ It , 

for i = L,R. Using (|57|l we may write the operator cor- 
responding to U 3 X : 

U? x - Ka = M^M^P+ClC^N, 
+P-Q + C^C*N 2 + Q-C&fi&Ns] , 
where we have defined the operators 



C' 



and 




Ml-- 

for i — i, R. The normalization factors are 

cj = 1 



(64) 



(65) 



(66) 



dim( P , g ) 
dim(p+i,g) 

dim(p,g) 
dim(p-i,g+i) 

dim(p,g) 
dim(p,g-i) 



(67) 



The C uv are operators that recover the Clebsch-Gordan 
coefficients as defined in Tab. |n] These were calculated 
using the "pattern calculus" of Biedenharn and Louck 
[2^ (see Appendix [BJ . Similarly, the N u recover the 
normalization factors defined in (|67|l . The operator cor- 
responding to U 3 X may be found using the identity 



U, 



IL , R = {-l) QL+QR (Ul VL _ Vjt r (68) 



where 



Qi = T? + \Yi + \ 



(69) 



and — vi = (Ti, —T*, —Yi). The phase factor follows from 
the standard choice of phase factors between states of 3 
and 3* |2J]. The desired operator is then 



U 3 ' «-> 

VLVR 



(—1 ^Ct + CRl/t 



(70) 
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TABLE II: Clebsch-Gordan coefficients for R ® 3 for SU(3). 
The tabulated formulas give the values for {p+Ap, q+Aq, T+ 
AT, T Z +AT Z , Y+AY\p, q, T, T z , Y; 1, 0, AT, AT Z , AY). For- 
mulas may be obtained by combining the results from Table 
[I] and Table ITTT1 where i — L, R is put in as appropriate. 



efficient 


Ap 


Aq 


AT 


AT Z 


Ay 


Formula 


°11 


1 





1/2 


1/2 


1/3 


r» 

J n ( -n 


f~li(b) 
°11 


1 





-1/2 


1/2 


1/3 


-^12 c 21 


^ii(a) 
°12 


1 





1/2 


-1/2 


1/3 


^ll c 12 


r ,i(b) 
°12 


1 





-1/2 


-1/2 


1/3 


ji i 

-I12C22 


ci 3 


1 











-2/3 


-<13 


°21 


-1 


1 


1/2 


1/2 


1/3 


- , 21 c ii 


s~ti(b) 
°21 


-1 


1 


-1/2 


1/2 


1/3 


-^22 c 21 


/-f*( a ) 
°22 


-1 


1 


1/2 


-1/2 


1/3 


Ti ~i 
J21C12 


°22 


-1 


1 


-1/2 


-1/2 


1/3 


ji i 
-<22 c 22 


C23 


-1 


1 








-2/3 


P 
1 2?, 


°31 





-1 


1/2 


1/2 


1/3 


-<31 c ll 


r<i(.b) 
°31 





-1 


-1/2 


1/2 


1/3 


ri i 
1 32 c 21 


(~ii{a) 
u 32 





-1 


1/2 


-1/2 


1/3 


Ti ~i 
J31C12 


°32 





-1 


-1/2 


-1/2 


1/3 


Ti n i 
-i32 c 22 


^33 





-1 








-2/3 


^33 



where the f refers to Hermitian conjugation with re- 
spect to the spin operators (e.g. (P ± ) 1 ' = P T , (T' i ± )t = 

T?, {Ciu) ] = &iu> ■ ■ ■ )• This result ma y be checked by 
directly constructing the operators in the same way as 
(|64|l with the proper 3* Clebsch-Gordan coefficients, and 
verifying that l|7U|l holds. 

The SU(3) Hamiltonian may be written 

H sv(3) =^i? 2 (r, M )-x Yl (Z(p) + Z\p)) 

r,/J. {plaquettes} 

(71) 

where 



z( P )= (-i) Si+S3 

Vl 1/2 V 3 Vi 

V„ lU2 (r, Li)V U2V3 (r + fi, v)Vl iVi (r + u, p)Vj 4 
E 2 \p, q) las = \{p 2 + q 2 +pq + 3(p + q)) |p, q) leg . 



and 



This is the final result for SU(3). 

Let us again find the required number of qubits and the 
dependence of the number of operations for a lattice size 
with M sites. The total number of qubits is ~ 8MDd, 
due to the eight registers required for each link each with 
approximately D qubits. With a similar arrangement 
of qubits as Fig. |3| we have a linear dependence of the 
number of operations with Md. The analysis for the 
worst case is identical to the SU(2) case, and thus we 
obtain a worst-case number of operations increasing as 
ex M 2 d 2 D. 



IV. PREPARATION OF THE INITIAL STATE 

In this section we discuss how the qubits should be ini- 
tialized in order to extract observables with respect to the 
low-lying eigenstates of the Hamiltonians. In particular, 
we have in mind the method of extractingthe energies of 
the Hamiltonians as discussed in Refs. [ill IT^ ] . In both 
of these methods, the qubits are prepared in an initial 
state with a non-zero overlap to the state of interest, as- 
sumed to be the ground state in this case. In the method 
given by Abrams and Lloyd, a sequence of operations in 
the spirit of phase estimation is performed on the qubits 
to extract the eigenvalues of the time-evolution operator 
U(t). In the method given by Somma et al, the expecta- 
tion value of U(t) is taken with respect to the prepared 
state, then a Fourier transform of measurement is per- 
formed on a classical computer to extract eigenvalues. 
The method of operating the time-evolution operator was 
discussed at length in Sec. IIHI However, it is clear that 
the initial state preparation must also be done efficiently. 

One potential choice of initial state is the strong cou- 
pling ground state, i.e. E(r,p;)\0) = for all r and [i. 
This is a good approximation for small values of x, how- 
ever, for large x the overlap with the true ground state 
becomes increasingly small. Since it is the weak coupling 
limit (x — > 00) that is generally of interest, it is clear 
that a better choice is necessary. In the following, we 
shall consider how this is done for the U(l), SU(2), and 
SU(3) lattice gauge theories. 



A. Initialization for U(l) Lattice Gauge Theory 

A popular choice of trial state for U(l) lattice gauge 
theory takes the form [2^] 



*io) = n 

pG{plaquettes} 



1 



exp [C{x)(Z{p) + Z\p))} |0) 

(72) 

where the plaquette operator Z(p) was defined in 114fl . 
|0) is the state defined by (|17|1 . and Af is a suitable nor- 
malization factor. C(x) is a variational parameter to 
be optimized with respect to the coupling x. This may 
be performed straightforwardly according to the method 
given in Ref. |2(|. 

Since the exponential operator appearing in (|72|l is 
non-unitary, one cannot simply write down the operator 
that prepares this state from the strong coupling ground 
state |0). It is still possible to prepare the state given 
in l|72|) by the following procedure. Concentrating on a 
single plaquette in the lattice, expanding the exponential 
we obtain 



*(p) 



Co + C 1 (z(p) + z\ P )) 



UZ{pf + z\p) 2 ) + 



(73) 
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For the U(l) case the expansion coefficients may be easily 
evaluated: 



o = E 



(C(x)y 



Un-l) 



Ii(2C(x)), (74) 



where // denotes the modified Bessel function of the first 
kind. The coefficients Q are kept until most of the dom- 
inant terms are found. For C{x) = 1, we find that the 
most of the dominant terms are contained in the range 
< I < Z max fa 5. 

Now let us see how we may prepare the state <|73|) using 
the formulation given in Sec. IIII Al There, we defined 
operators L ± that increase and decrease the \l) register 
by ±1. Let us introduce a similar operator acting on a 
register with the following properties: 



L(I') 



10, 
o 



in = o 

otherwise 



(75) 



The Hcrmitian conjugate operator to this is 
L\l')\lU- 



|0) r 




if / — /' 
1 . ■ (76) 
otherwise v ' 



These operators are restricted versions of the operators 
L^, that are specialized to shifting between the states 
|0) and |Z). An example of these operators is given in 
Appendix [X] We may construct the restricted versions 
of the plaquette operator as follows 



Z«(p) = L(l) L(l) L\l) L\t), 



(77) 



where we have labeled the links that each of the opera- 
tors apply underneath each operator. Now consider the 
operator 



U, 



z( P ) _ 



exp 



0(Z,(p)-Z/(p)) 



which is a unitary operator. The advantage of working 
with these restricted versions of operators is that the re- 
sult of acting these operators may be written in closed 
form: 



cos (f>\0) + sin 
Uf (p) |Cr") = { cos <£!□'> - sin0|O) 



if V = 
if V = I , 
otherwise 

(79) 

which is easily verified by expanding with respect to 4>. 
In (|79)l we defined 



Zl(p)\0)r 



(80) 



The initial state preparation is then performed by the 
operation 

n flV^o^oMiow (si) 

pG{plaquettes} l — \ 



The values of the <j>i are fixed by the Cz calculated previ- 
ously. For example, on a single plaquette with i maJI = 1, 

£/f (P) (MU^ P) (4>-l)\0) = COS 0! COS 0_!|O) 

+ sm0_ip- 1 ) +sin0icos0_ip 1 ), (82) 

and thus cos tpi cos <f>—i = (q/\/J7, sin</>_i = 
sin 0i cos0_i = Cx/vJ\. This is then repeated for every 
plaquette in the lattice, which completes the initializa- 
tion procedure. 

B. Initialization for SU(2) Lattice Gauge Theory 

A suitable trial wavefunction for SU(2) lattice gauge 
theory is H3 



*|0> 



n 



i 



p£{plaqucttcs} 



exp[C(x)Z(p)}\0). (83) 



Concentrating on a particular plaquette, we may again 
expand the exponential to obtain 



*(p) = -= Co + CiZ 1/2 (p) + c 2 Zi(p) 



(84) 



where 



Zj(p) = Tr \u j (r, fi)U j (r + n,v)U j \r + v, n)U j \r,v) 

(85) 

and j labels the representation of the U operator. The 
terms in the expansion (|84|l may be calculated with the 
help of the identity 

Z 1/2 ( P )Z j (p) = Z,-_ 1/2 (p) + Z j+1/2 (p). (86) 

In analogy to the U(l) case, let us define restricted ver- 
sions of the and ML operators as follows: 



J{j')\j, m L ,m R ) 
M L (m')\j, m L ,m R ) 
M R (m')\j,m L ,m R ) 



\j',m L ,m R ) IOS if j = 
otherwise 

\j, m', m R ) Teg if m L = 
otherwise ' 

\j,m L ,m') Teg if m R = 
otherwise 



The restricted plaquette operator is then defined as 

Zjm 1 m2m3irii(p) 

J(j)M L ( mi )M R (m 2 ) J(j)M L (m 2 )M R (m 3 ) 

s ^ ' 

x J(j)M L (m 3 )M R (m 4 ) J{j)M L {m A )M R (m 1 ){^7) 



where we have again labeled the links that each operator 
acts on underneath each operator. We now define the 
unitary operator 



tjZ(p) 



(P)" 



m 

(88) 
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which has the operation 



U 



zip) 



\n 4 



COS 0| 0) 

cos d>\ □ 



sm < 



P 



sin 0| 0) 



if f = m' t = {) 

if f = j, m'i = 
otherwise 



(89) 



where we defined 

□ 



jmim2fn3'm4 / 



J jmim2m3rti4 



(p)|Q>- 



(90) 



The initial state is then prepared by a product of these 
operators on the strong coupling ground state: 



17 17 17 u z(p) u- 

11 11 11 LJ 3mim2m 3 m 4 lr3mim2m3m4 

p£ {plaqucttcs} j — 1/2 mim^m^m 4 



)|0) I 



where each of the m% summations run from — j max to j„ 



(91) 



The 



are fixed using the Q coefficients in a 



similar way to the initialization for U(l) lattice gauge 
theory. This initialization is then carried out on each 
plaquette of the lattice. 



operations required for each JJ z( ^p) (<^) by a factor of the 
number of qubits used, which is equal to DdM . At worst, 
we therefore find that the total number of operations is 
proportional to Dd 2 M 2 . 



E. Other choice of initialization 

Our analysis in the previous subsections were restricted 
to considering particular forms of trial wavefunctions. 
These are clearly not the only trial wavefunctions that 
may be used, for example, correlated plaquette terms 
may also be included p9l | . These trial wavefunctions may 
also be constructed using the methods given in this sec- 
tion. Inclusion of states that have long-distance correla- 
tions over the lattice will naturally increase the complex- 
ity of the initialization procedure by a factor proportional 
to the size (e.g. the length of the perimeter of the Wilson 
loop operator) of the state being created. Since this is 
bounded by the size of the lattice, the total number of 
operations will still be ~ poly (M). 



MEASUREMENTS 



C. Initialization for SU(3) Lattice Gauge Theory 

The initialization procedure for the SU(3) case is very 
similar to the U(l) and SU(2) cases. We will therefore 
only show a sketch of the results. A suitable trial wave- 
function for SU(3) is [H 

*|0)= -^exp[C( iB )(Z(p) + Zt(p))]|0). 

pG{plaquettes} ' 

(92) 

In analogy to the U(l) and SU(2) cases, we de- 
fine restricted plaquette operators ^AiA 2 a 3 a 4 , 
which creates the state associated with the term 
U \ 1 \ 2 U \ 2 \ 3 U \3\ i U \!\ 1 in the plaquette operator Z(p). 
The initialization is then performed by a product of 
unitary operators defined using the restricted plaquette 
operators: 

II II ^ilA3A 4 (^A 1 A 2 A3A 4 )|0) rcE . (93) 

p£{plaquettes} 0jAlA 2 A3A4 

D. Efficiency of the initialization procedure 

Let us now estimate the dependence of the number of 
operations of the initialization procedure with the lattice 
size M. The total number of times the U z ^{<j)) operator 
must be applied is clearly proportional to the number of 
plaquettes in the lattice oc dM. Therefore, with a layout 
of qubits that preserves the local nature of the operation, 
the number of operations will be oc Md. However, an 
unfavorable layout of qubits can increase the number of 



To obtain results from the simulation, one must be 
able to extract quantities such as the expectation value 
of operators and the spectrum of the Hamiltonian oper- 
ator. The energy spectrum of the Hamiltonian may be 
obtained following the method outlined in the beginning 
of Sec. IIVI The most interesting expectation values to 
be evaluated are of the Wilson loop operators, defined as 

W r = Tr J] U{r,n), (94) 

where T denotes the contour in the lattice that the Wilson 
loop follows, and U(r, /i) denotes the unitary operator 
defined for the link labeled by (r, /x) for SU(iV). We may 
immediately transcribe these operators into those using 
spin operators using the results of Sec. IIIII For example, 
for SU(2), 

W T = ^2 Vm im2 (ri,Hi)V m2m:3 (r2,fJ.2) ■ ■ ■ 

mi 2 , . ..,mjv 

(r N ,/j, N ). (95) 

Analogous results follow for the U(l) and SU(3) cases 
hold. Expectation values of this operator may be found 
for example by the general phase estimation method 0]. 
This involves performing the operation 

U(t)=e~ lWr , (96) 

which may be performed in a similar way to the methods 
explained in Sec. IIIII 
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VI. SUMMARY AND CONCLUSIONS 

We have reformulated the U(l), SU(2), and SU(3) lat- 
tice gauge theories into a form that can be simulated 
using spin qubit manipulations on a universal quantum 
computer. This was done by constructing spin operator 
versions of the link operators E and U appearing in the 
Hamiltonian versions of the theories. Each link on the 
lattice is represented by a set of registers that keep track 
of the state of the link. In all cutoff corresponding 

to the maximum representation of the group was found to 
be necessary to keep the total number of qubits required 
for the simulation finite. This is necessary as even a sin- 
gle link has an infinite dimensional Hilbert space. These 
high representation states have an energy that grows as 
the square of the index that labels the representation (e.g. 
j(j + 1) for SU(2)), and thus the impact of truncating 
these states can be carefully controlled, as has been found 
in past studies. For the SU(2) and SU(3) calcu- 
lation of Clebsch-Gordan coefficients is required for every 
application of the plaquette operator Z(p). Explicit for- 
mulas were given for the Clebsch-Gordan coefficients for 
SU(2) and SU(3), which are precalculated on a classical 
computer for a given j max or p max , g max respectively. This 
is then written in operator form which requires an over- 
head corresponding to ~ poly(j max ) or ~ poly(p max , g max ) 
respectively during the simulation. The number of qubits 
necessary to keep track of a lattice with M lattice sites 
and dimension d was found to be ~ dDM, where D is 
the number of qubits required to keep track of a single 
link state. The number of operations required to perform 
the Hamiltonian evolution of the system was found to be 
proportional to between a linear and quadratic function 
of the total number of lattice sites. This is dependent on 
the architecture of the quantum computer that is avail- 
able. The simulation may therefore be said to be effi- 
ciently implementable in this sense. Other aspects of the 
simulation, such as the initialization of the qubits into a 
suitable configuration for extracting low-energy proper- 
ties of the system were also considered and were found to 
have a similar dependence with lattice size, and therefore 
does not spoil the efficiency of the whole procedure. 

With the current technology available, the kind of sim- 
ulation that we suggest in this paper is clearly still a long 
way off. For example, a typical size of lattice for a QCD 
simulation today has the order of 10 4 lattice sites. By 
our scheme for SU(3), we would require ~ 10 5 D qubits, 
where D is the size of the register dependent on the cut- 
off that is employed. Note that this does not include any 
quantum error correction, and therefore the real figure 
will be larger than this. However, the method here has 
the advantage that it is virtually an exact calculation, 
up to the cutoff on each link which is controllable such 
that its impact is negligible. All the states associated 
with the low-lying excitations are included. With the in- 
clusion of the fermion fields, which adds a further ~ M 
qubits to the total, this would allow a direct and exact 
simulation of QCD itself - an unthinkable prospect on a 



classical computer. Although significant advances have 
been made in the lattice QCD programme in the past 
few years 0, such an "exact" simulation may be a use- 
ful check on the results of the classical simulations and 
may even offer new insights into the physics of QCD and 
related models. 

We have restricted our focus in this paper to pure 
gauge Hamiltonians, i.e. no fermions. For a realistic 
simulation of QCD, we also require fermions in addition 
to the gauge degrees of freedom that were examined in 
this paper. One way that fermions may be included into 
the formulation is to transform the fermion operators to 
spin operators via a Jordan- Wigner transformation 0. 
This involves long strings of operators that are of order 
the size of the lattice M. The space required for includ- 
ing the fermion adds M more qubits to the entire lattice, 
which far less than the number required for the gauge 
fields. This should therefore not affect the main results 
of this work. 

We note that the same methods as that given here may 
be straightforwardly generalized to SU(iV) lattice gauge 
theories, and we expect that these may also be efficiently 
implemented. One vital ingredient that makes this pos- 
sible is due to the Clebsch-Gordan coefficients for the 
product of the fundamental representation of the group 
to an arbitrary representation of the group according to 
the "pattern calculus" of Biedenharn and Louck . 

There are many ways that the implementation that we 
describe in this paper may be improved, particularly with 
respect to the total number of qubits. In particular, there 
are many wasted states on a particular link associated 
with states that do not exist. For example, for the SU(2) 
case, unphysical states such as \j,rriL > j,mR > j')r eg 
are included, as far as the qubit implementation is con- 
cerned. Likewise, unphysical gauge variant states are also 
included. We chose to ignore such inefficiencies as this 
would introduce further complexities that are unneces- 
sary for the main goals of this paper, although these 
should be possible, if necessary. 
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APPENDIX A: SOME SIMPLE REGISTER 
OPERATORS 

To illustrate the methods used in Sec. IIIII let us con- 
sider a simple example of the operators that could be 
used to manipulate the registers. Considering the U(l) 
Hamiltonian Ij26(l , we see that this consists of operators 
L ± and E 2 that are defined according to (|2"2")) - l(2"5|) . A 
simple, but qubit inefficient way that this could be im- 
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plemented is as follows 

imax-l 

;=-; max 

L = E (A2) 

i=-Z max 

E ;2 K + !)A (A3) 

i=-Im,x 

where the and a z are Pauli spin operators. These 
operators are designed to act on states with a single up- 
spin: | H . . . f . . . I). The location of the up-spin keeps 
track of the register: 

\l) res = <rt\il...i). (A4) 

A similar type of register was considered in Ref . 0] . Sim- 
ilar operators are used in the SU(2) and SU(3) Hamilto- 
nians. 

The operators that recover the Clebsch-Gordan coef- 
ficients may be constructed by a linear combination of 
diagonal operators that act on the registers. Considering 
the SU(2) case, let us construct the operator associated 
with the coefficients c mn . We may write 

Cmn = ao + a x J z + a 2 M z + a 3 (J z ) 2 + ..., (A5) 

where 

J Z \j)r. S =j\j)r. S 

M z |m) rog = rn|m) reg 

and the a, are coefficients to be determined. The number 
of terms in (|A5|) is determined by the number of distinct 

I 



/ miN m 2N ■ ■ ■ 

miN-i 

m 12 

V "in 



The integers rriij satisfy the "betweenness" conditions 

my+i > rriij > mj+y+i (B3) 

The corresponding state in the \J, M) or \p, <?, T, T 2 , Y) 
notation for a particular Gelfand pattern may be found 
by the following correspondences - for SU(2) [23j 

J = |(mi2 -m 22 ), (B4) 
M = mn — |(mi2 +77122), (B5) 



(j, m) combinations. For a j max = 1, there are three 
(j, m ) states, and thus there are three terms in (|A5|) . In 
this case, the coefficients for en will be ao = 1, a% = 
1 + 1/V2, and 02 = 1 - 1/V2. 

Restricted versions of the operators considered in Sec. 
IIVI (i.e. Eqns. (|73|l and l(7B)l) may be constructed as 
follows: 

L(l) = <r a+ (A6) 
LHl) = <j+ a-. (A7) 



APPENDIX B: CALCULATION OF 
CLEBSCH-GORDAN COEFFICIENT FORMULAS 

Clebsch-Gordan coefficients in Tables HI - IIIII are calcu- 
lated using the "pattern calculus" method developed by 
Biedenharn and Louck [2^, . The method allows the 
calculation of Clebsch-Gordan coefficients for the combi- 
nation of an arbitrary state of SU(iV) with the fundamen- 
tal representation of the group. For example, Clebsch- 
Gordan coefficients for j ® \ in SU(2), or R® 3 in SU(3) 
may be found. 

Each finite dimensional irreducible representation is 
specified by a set of N ordered integers (positive, zero, 
or negative): 

[m] N = [m 1N m 2N . . .m NN ] (Bl) 

where m\N > nri2N > • • ■ > tunn- A state in this ir- 
reducible representation is represented by a triangular 
array called a Gelfand pattern: 



(B2) 



For SU(3) we use the relations 



p = 


77113 - 77723, 


(B6) 


q = 


77723 - "733, 


(B7) 


T = 


i(777 l 2 - 777 22 ), 


(B8) 




77111 - 5(77712 + 77122), 


(B9) 


Y = 


77112 + 777 22 - | (mi3 + 777 2 3 " 


K77733). (BIO) 



Using the above relations, the method given in Sec. Ill 



mjV-lJV 777ATJV \ 

TOJV-UV-I 



777 22 
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of Ref. |22( may be used to generate the matrix ele- 
ments of reduced Wigner operators. For SU(2), matrix 
elements of the reduced Wigner operators are equal to 
the Clebsch-Gordan coefficients for SU(2). For SU(3), 
the matrix elements of the reduced Wigner operators are 



the isoscalar factors for SU(3). The formulas can be cal- 
culated up to a ± sign. The sign is fixed by comparison 
to tables of Clebsch-Gordan and isoscalar factors (see for 
example Ref. pjof). 
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TABLE III: Isoscalar factors for i?(g)3. Isoscalar formulas are 
obtained by substituting ttf = 4p + 2q ± 6T, - 3V"» + 9 ± 3, 
nf = 2p-2q±6T l +3Y i -3±3, nf = 2p+4q±6T i +3y i +3±3, 

r 1 = (i+p)(2+ P + g ), r 2 = (i+p)(i+ g ), r 3 = (i+ q )(2+ P + q ), 

Ti = 432(1+Ti), T 2 = 432Ti, and T 3 = 36 into the tabulated 
formulas. The label i — L, R is put in as appropriate. The 
reduced Wigner operator used to calculate the formulas are 
also shown. 



Symbol Reduced Wigner operator Formula 
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